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Materials with negative permittivity and permeability can overcome the diffraction limit, thereby making 
the sub -wavelength imaging possible. In this study, we analyze the effects of gradient index on a half-infinite 
perfect lens. We assume that the sharp interface between the vacuum and the negative-index material is 
replaced by a smooth transition profile such that the index gradually changing from positive to negative. 
Interestingly, we find that if the graded index profile is modeled by a tanh function, we can have closed-form 
analytical solutions for this problem, which is a distinct advantage as numerical solutions are not accurate 
for evanescent waves with large transverse wave vectors. By analyzing the analytical formulas we confirm 
that a nonzero total absorption can occur even for a near-zero absorption coefficient in the steady-state limit 
and the image plane contains multiple sub -wavelength images of an object. 

Pendry 1 proposed that a slab with e = n = — 1 can overcome the diffraction limit leading to the perfect lens 
concept. The mechanism is that the perfect lens can amplify the evanescent waves and restore the high- 
spatial-frequency information of an object. The effects of absorption on perfect imaging of a finite-length 
slab were adequately considered 2 ' 3 . Several experimental demonstrations on the subwavelength imaging system 4,5 
have been achieved. 

In practice, a negative-index material cannot have an infinitely sharp interface between the material and 
vacuum. From an effective medium point of view, the transition from vacuum to negative index should be a 
smooth function. In that sense, the perfect lenses can be modeled as positive-to-negative transition materials 6 . As 
such, they have been modeled as layered optical materials with gradient optical indexes e and/or [i continuously 
changing from positive values to negative ones. In Ref. 7 and 8, it is found that near the transition point for £, the 
fields for oblique incidence present large enhancement with enhanced absorption. The absorption is found to be 
nonzero even in the lossless limit. In Ref. 6, similar effects are obtained when both t and /.i linearly pass the same 
transition point. The question we are addressing here is the effect of the zero-index point on the quality of perfect 
imaging. 

In order to answer that question we consider a simplified case, a half infinite negative-index material with a 
gradient index profile described by the function tanh. For simplicity, that material is free of frequency dispersion. 
To study the imaging effects, we need to analyze the whole spectrum of the object. However, numerical methods 
(such as transfer matrix method) for evanescent waves can be challenging as the accuracy depends on the 
cancelation of large matrix elements (as we discuss later). Fortunately, it turns out that the tanh function model 
enables us to analytically treat the transmission and reflection of propagating waves as well as properties of 
evanescent waves. Therefore, the analytical results can provide us with a precise description of imaging and 
absorption in these novel systems. 

Before we systematically investigate the imaging properties of this structure, we confirm the existence of the 
large absorption near the zero-index point according to Ref. 6. The maximum absorption and several asymptotic 
properties are also addressed. We then analytically and numerically answer the question posted above about the 
effects of large absorption (arising from the zero-index point) on the perfect imaging. The main feature is that the 
image plane contains multiple subwavelength images of an object and the distance between the subimages is 
proportional to the characteristic transition length. This is a unique feature attributed to the gradient index of the 
perfect lens. 

We note that this study provides an example of transition materials" -13 that can be treated fully analytically even 
for oblique incidence. 
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Results 

Material model and its analytical solutions. The permittivity and 
permeability we use to model a gradient- index passing from positive 
to negative index are e — (l= tanh(px + id). We note that £ = /(—> — 1 
for x — » — 00 and e = fl— > + l for x — » +°o. The parameter p 
characterizes the gradient of the transition. The imaginary part of £ 
and p is positive for 0 < 5 < nil, and therefore the material is a 
passive absorbing material under that condition. One example is 
shown in Figure 1. For simplicity, we set the velocity of light c = 1 
and the vacuum parameters £o = 1 and p 0 = 1 in the whole study. As 
we treat all quantities as dimensionless, the results can be applied to 
different length scales after some rescaling. 

The Helmholtz equation for TE waves have already been written as 
in [6], 

d 2 E 1 dp dE 
dx 2 p dx dx 



I 2 2 

[erpf - 



b 2 )E = Q. 



(1) 



where a = oj/p, b = oj sinO/p with a) and 0 being the frequency and 
angle of the incident wave. This equation can be transformed into 
hypergeometric equation after changes of variables, y = p 2 and E = 



(1 — y) A F (y) where X ■ 

d 2 F 
~df~ 



- b 2 . The transformed equation is 



F = 0, 



2X+1 dF 



l-y dy 47(1-7) 



which is a hypergeometric equation with indexes a = X + ia/2, f} = X 
— ia/2 and y = 0. 

Analytical solutions for this material. In the following Methods 
section, we show all required steps to obtain the final analytical 
results. These steps generally involve the properties of hypergeo- 
metric equations and series. Here we only refer to the final results 
which are summarized in Eq.24. 

Basic properties of this material. In this part, we study some basic 
properties of the gradient- index material by investigating the 
analytical solutions. It includes the discussions about maximum 
total absorption rate and asymptotic behaviors for some limiting 
cases. These properties are helpful for us to understand the effects 
on imaging. 

Maximum total absorption. As discovered in Ref. 6, the total absorp- 
tion of this positive-to-negative transition material is nonzero even 
when the absorption coefficient is near zero. In this section, we 
confirm and develop this result for our specific model. 

Before we investigate the total absorption in the lossless limit 5 —> 
0, we first prove the identity |7c| + |T| = 1, or equivalently 
I C (1) I = 1 + I C {2) I . We can easily see that | C (1 _> | = C ( _l> > 1, since oc 
and p are purely imaginary. Then we calculate I C': 2 ' I . 



I c(_ 2) f=cw c<_ 2) * = cL 2) («,/*) c(_ 2) ( - 



-/?) = 



4sin 2 7t/isin 2 7ta 
sin 2 Tt(a + /?) 



(3) 
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Figure 1 | The permittivity c of the material, p = 1, 5 = 0.01. 



Then follows the identity | | = 1 + | | . 

Given the above identity, the total absorption A = 1 — \1Z\ 2 — |T| 2 
reaches its maximum 1/2 under the condition \1Z\ = \T\ = T = 1/2. 
Then we finally arrive at the incident angle for maximum absorption^ 
which is given by 

In cosh 710J / p 



cos 0 = ■ 



TZOj/p 



(4) 



That is to say, even in the lossless limit <5 — » 0, for any given fre- 
quency, this material can present maximum total absorption 1/2 at a 
specified angle given by the above equation. 

1Z and T in several limiting cases. We only consider the limit <5 — » 0 in 
this section and do not introduce frequency dispersion to our model. 

For oblique incidence 0^0 and high frequency CO—* + °°, it can be 
easily proved that j C'j 2 ' j -> + 00. Therefore lim (0 _ ( + . x |T| =0, and 
limoj^ + oolT?.! = 1. That is, for high frequencies, the light is totally 
reflected. However, this is the case when no dispersion is taken into 
consideration. 

For oblique 0 7^ 0 but low frequency oj — > 0, we have a. 0 and /i 
-h> 0, foUowed by | C^> | = 0, |cW| = l and lim^ol^l = 1, 
lim m ^ 0 |7?-| = 0- F° r l° w frequencies, the material behaves like a com- 
plementary material with total transmission. 

For normal 0 = 0 and a finite to, we can easily get limg_, 0 |'?"| = 1> 
\img^ 0 \TZ\=Q. That means for normal incidence, we have total 
transmission which is the characteristic property of complementary 
materials suggested by Pendry and co-workers. 

For large incidence angle near n/2 and a nonzero cu, we can prove 
that Ic' 2 ']— > + co. Then follows ]img_ >7t /2|'7 _ | = 0, and limg^ n / 2 |^| 
= 1. For glancing incidence, the light is totally reflected. 

We should remind ourselves that for a sharp interface between the 
vacuum and negative-index material without smooth transition, the 
light will always be totally transmitted. We can see that for the gra- 
dient index, not only can we have a transmittance less than unity, but 
also have a nonzero reflectivity. Besides, the total absorption is non- 
zero even in the lossless limit <5 — » 0 as discussed in the previous 
section. 

Comparison with numerical calculations. In order to verify the 
analytical formula, we numerically compute the total absorption 

Total absorption for<o=27i, p=1, 8=0.01 



numerical results 
analytical results 




Figure 2 | The comparison of the analytical and numerical results for the 
total absorption under m = In, p=\, and 3 = 0.0 1 . The numerical results 
are obtained by use of transformation matrix technique. The material is 
confined within [ — 5, +5], and is discretized into 3000 uniform layers. 
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(l-\Tl\ -\T\ ) of this material by use of the transfer matrix 
technique. The analytical and numerical results are compared in 
Figure 2 under the condition that w = 2n, p = 1, and 3 = 0.01. 
We can see that the two different methods give the same results. The 
total absorption can exceed the 1/2 bound because in this case we 
consider a finite absorption parameter S = 0.01. 

For evanescent waves, the numerical method based on the stand- 
ard transfer matrix is inaccurate, because the nonzero transmission, 
which is necessary for perfect imaging, relies on the accurate cancel- 
lation of very large matrix elements. Therefore, if we want to accur- 
ately study the whole spectrum, analytical is the best tool. 

Effects on perfect imaging. In Ref. 1, Pendry proposed that a slab 
(perfect lens) with e = p = — 1 can overcome the diffraction limit. 
This is made possible by the fact that the perfect lens can ideally 
restore the information contained in the evanescent waves. 
However, for gradient-index material, the reflection and 
transmission can be very different from the sharp-interface case as 
we showed in previous sections. In the following, we analyze the 
effects of a graded transition from positive to negative index. 

In the sharp-interface limit, p —* +°°, the transition materials 
approaches a half-infinite perfect-lens slab. A perfect image of an 
object in another side of the interface is therefore possible since the 
transmission for all transverse wave vectors are unity. 

However, as can be seen from previous sections, if the transition is 
smooth (as in Fig. 3), i.e., p is finite, the transmission is not unity even 
for the propagating wave. If we set a monochromatic line object (not 
a point) in the plane x = — a, with an arbitrary transverse amplitude 
profile E{x = —a;z) and Fourier components E(x= — a;k) = 
E(x= — a; z)e' kz dz, we will obtain the field amplitudes in the 



image plane x = a, 

E{x = a;z) 



2n 



E(x= -a;k)T{k)e-' kz dk 



where, according to the results in Methods section, 

sin 2*; 



-cosh Tzwj p ' 



(5) 



(6) 



with a = — ^/k 2 — to 2 J 2p being the principal value of the square 

root. We have also assumed that 8 — » 0. 

Because T is a periodic function of a, we can obtain the discrete 
Fourier expansion of T, 

+ x 



A remarkable feature of each Fourier term exp (i%n^ (k 2 -co 2 )/ 
p) is that for large k, it approaches a periodic function and thus 
contributes a resonance when inversely Fourier-transformed. Let 
us approximate exp inn^J (k 2 - 



o 2 )/p~exp im Wf> for large k. 
Then the images in the real z axis can easily be obtained after inverse 
Fourier transformation. 



n — -t ou i. + co 

E(x = a:z)~ a n \ E{x= - a;k)e- ikz+i \ k \™> <f <dk. (8) 

n — — ao — co 

If we consider a simple Lorentzian object E(x = —a;z) = ~, 

nz 2 + A 

the image should be 



^ a„ nn p-iA 
E(x = a;z)~ > 



nfr'oo 71 Z 2 -(nK/p — z'A) 



(9) 



T-- 



exp 



(7) 



The above result indicates that the image plane contains multiple 
peaks. The widths of the all peaks remain the same as that of 
the object which is a Lorentz spot. The distance tween subimages 
is nip which means that the emergence of multiple images can be 
attributed to a finite p or the gradient index. The amplitude a n of 
high order images exponentially decays with a characteristic order 
1 

fir — . 

In cosh(7rw/p) 

In order to obtain more accurate results for the image profile, we 
numerically integrate Eq.5 for the parameters <5 = 0, A = 0.02, p = 5, 
co — 1. These parameters means that the characteristic width of the 
transition layer is about 1/30 of the wavelength. The width of the line 
object is about 1/300 of the wavelength. The numerical results are 
displayed in Figure 4. That is to say, the width of the line object is 
much less than that transition width of the layer which itself is still 
much less than the wavelength. We can see that image profile appears 
in such a way that many peak emerges just as sub images of the line 
object. We can also verify that the distances, about n/5, between these 
peaks can be predicted by the Eq. 9 quite well. The widths of these sub 
images are comparable to that of the line object, which is also in 
accordance to the Eq. 9. 

Let us study a more complex object, E(x = —a;z) = 
A 1 A 1 

s r H s r , which contains double 

n( z + d/2) 2 + A 2 n( z -d/2) 2 + A 2 

Lorentzian objects. Following the same procedure, we obtain the 
corresponding image 



with a„ — Ux 



■1)X "for« = l,2,- ■ -,a 0 = \% 2 \a_ x = \x 



0 for n= — 2, — 3, • • •. Here x = cosh(7rco/p). 




-5 -5 

Figure 3 | The setting of the material. 




-10 -5 0 5 10 

z: Transverse position in the object (image) plane 

Figure 4 | Single line. The field amplitude distribution in the object plane 
(red) and image plane (blue). The relative strength between the red and 
blue curve does not have meaning. 
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u n 



nn/p — iA 

if- 

nn/p — iA 



K(z-d/2) 2 -(nn/p-iA) 2 



(10) 



(z + dj2) l -{nn/p-iAy 



We can see that two sets of multiple images overlap with each other. 
We can confirm the approximation in Eq.(10) by direct numerical 
calculation, as shown in Fig. 5. 

Discussion 

We have not included the finite-length effect (finite-length slab) in 
this study, and thus can not address the frequency cut off proposed in 
Ref. 3. If finite length, gradient index and absorption are all consid- 
ered, we can expect that we still see multiple images but the detail 
profile of each subimage needs to be investigated. Since the material 
we discuss in this study is free of dispersion, we cannot address the 
formation of the image in the time domain. To analyze the more 
realistic material with dispersion we may consult to numerical meth- 
ods and full analytical results are not possible. However, we can 
expect that it may take a relatively long time to form the image 
because of the field enhancement near the zero-index. These are to 
be studied in the future. 

In summary, the transmission and reflection properties of a 
graded positive-to-negative materials can be very different from 
the ones with idealized sharp interface. These differences are mostly 
attributed to the large field enhancement near the zero-index point. 
To study the effects of graded index on perfect imaging, analytical 
solutions are obtained for smooth transition material with a "tanh" 
profile. The analytical expressions provide us insights into this prob- 
lem. Expressing the fields in terms of hypergeometric functions, we 
obtain the closed-form analytical expressions for the reflection and 
transmission coefficients. The results are confirmed by direct numer- 
ical method. Our analytical results uncover the smooth transition 
effects on imaging: emerging of multiple images. The distance of 
these subimages is shown to be directly related to the gradient of 
the transition. These effects can be important for understanding the 
subtle properties of the perfect lens. 

Methods 

Solutions of the Helmholtz equation for x — » +oo, x — * — 00 and — oo < x < +oo. 
According to the material model we have adopted, for x — » ± oo, we have fi — > ±1 and 
y — » 1. As shown in the following section about Kummer solutions, the field can be 
expressed as linear superposition of two basic solutions, 



E ± (x) = cf (l-ffa + dg (1- M 2 )V, 



(11) 



£ 0.8 



T3 
CD 
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-10 -5 0 5 10 

z: Transverse position in the object (image) plane 

Figure 5 | Image profile for a source of double line. The field amplitude 
distribution in the object plane (red) and image plane (blue). The relative 
strength between the red and blue curve does not have meaning. We use d 
= 10A, and other parameters are the same as in Fig. 4. 



where w 3 4 , as well as w 12 below, are Kummer solutions for hypergeometric equations 
(shown below). 

For — oo < x < y is near the zero pole. The two basic solutions are, 

{\-fff and (1-,, 2 )V (12) 

where w 2 , defined in Eq. 26, is one of the Kummer solutions near the zero pole and/is 
the solution with logarithm since we are working with a degenerate hypergeometric 
equation with index y — 0. 

Similar to the expression of Neumann functions in terms of Bessel function, the 
logarithm solution/can be expressed as limit of linear superposition of Wi and w 2 of 
nonzero y index, 

Wi 



/=lim 



r(y) 



-a/?W2 



1 



(13) 



which will be proved in the Methods section. Specifically, according to the L'Hospital 
rule, 



/= lim- 



r(y) 2 



-otj3w' 2 



r'(y) 
r(y) 2 



(14) 



dwi dw 2 d 

where w i — — — , w 2 — —— — — w 2 In y-\-y ' — 2 Fi(ot — y + l,p — y + 1; 2 — y\y) 
dy dy ay 

with 2 F X being the hypergeometric function. 

However, if we use the principal values of f(y), the field is discontinuous 14 along the 
real x axis at zero point x — 0, corresponding to [X — tanh(/<5) — i tanh 5. Therefore, for 
0 < x < +=c and — °° < x < 0, we should use different superpositions of w 2 and /to 
make sure that the field is continuous. 

The discontinuity of the principal value off isf(x~ 0+) — f{x— 0 — ) — 2y.[iniw 2 {x 
= 0) according to Eq.(14). Then, if the field for 0 < x < +cc is 

E + (x)^A(l-fi 2 Yf + B(l-fi 2 Yw 2 , (15) 
the field in — °° < x < 0 should be 

£_(x)-A(l-/( 2 )V+(B + 2a/57rfA)(l-/( 2 ) A w 2 , (16) 

to ensure the continuity across x — 0. 

Relations amongC^ 2 ^ and A, B. The different field expressions {11, 15,16) should be 
equal in equivalent regions, which imposes connection relations among the coeffi- 

- 0 in x + oo region, according to the Kummer connection 



cients C ( + 2) and A, B. 

If we set d^UCf-- 
formula between w 3 and w 1)2 in Appendix A, as well as Eq. (13,15), we have 

r(y-l)r(a+J?-y+l) 



r(l-y)/>+0-y+l) 

W3 — — , , — W\ - 



- H'2 



(17) 



/>->.+ i)rQ3- y +i) " ' r(a)r(fi) 

= im -A(y)4ir(y)]w 2 +A(y)w 1 . 
with B — lim 1 ,_ (0 B(y) > A = ]im Y ^oA(y). After straightforward calculation, we have 
4 _r(a+/i+l) 

(18) 



a/Jr(a)r(/J) 
1 



i//(l) — 0.577215664.... and ijj is the derivative of the logarithm of V 



where y E 
function 

Then according to Eq.{15,16,17), in x — » — 00 

E_ (x) - (l -f-i 2 ) A w 3 +2a07TiA(l - n 2 ) A w 2 . 



(19) 



By use of Kummer connection formula between w 2 and w 3 4 in Appendix A, we have 



AajS(2ro)r(-a-$l AaP(2ni)r(x + P) H 2 ,>_ 



)F(l-« 



which, compared with Eq.(16), gives 



C« = l + 



(2ni)T(-x- 



r(i-a)r(i-« r(«)r(/5) 



r(a + /5 + l) ^ . 2 sin no. sin 7tji 



sin 7c(a + / 



p, = {2ni)T(a + P) r(ot + ^+l) 

- r(c+i)r(/;+i) r(«)r(/?) 



(20) 



(21) 



Reflection and transmission coefficients. Now we have obtained continuous 
solutions along the real x axis. For x — > E + — (1 — ^ 2 ) A w 3 ; for x — > — r - 
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£_ — C^) (1 — /i 2 ) A w 3 + (l — f.t 2 ) A w 4 . According to Appendix A, the asymptotic 
behaviors for x — > or^ ~ 1 are that tv 3 ~ 1 and w 4 ~ (1 — y) -2 ^. 
Therefore, for x 

E + ~(i-ff^pV^ e ^sV?=V 2 -iV^^ (22) 
which is an outgoing wave; for x — » — °°, 

e _^q{1) e -ixpjai - V e + <W<* 2 -b 2 2-Wa 2 + q{2) ^xp^/a 2 -b 2 ^^a 2 -b 2 2 Wa 2 -b 2 ^ ^3) 

where Ql^ term represents incident and reflected wave, respectively, since the 
energy flux and wave vector point in opposite directions in this double negative 
region. 

From Eq. (22,23) above, we can directly read out the reflection 71 and transmission 
T coefficients, 



11= 



C (2) 



- 0 - 2,3 vV - b 1 + i2 In 2 vV - b 1 



(24) 



' C<i> 



Kummer solutions for hypergeometric equations. For completeness, we include 
some information about hypergeometric equations used in this study. One can also 
refer to mathematics handbooks. The hypergeometric equation is 



, d 2 w r , , , dw 

z(l-z)-^ + [r-(l + a+j8)*]^-a0tv=O. 



(25) 



When none of y, y — a — [i, a — [5 is an integer, the two basic solutions near 0 < \z\ < 1 
are, 



w 2 {z) Vi (a- V + I.P-Y+ 1; 2-y; z); 
the two basic solutions for 0 < \z — 1 < 1 are, 
>v 3 (z)= 2 Fi(og5;a + ;8-y+l;l-z), 



w 4 (z) = (i-z) 7 a ViCv-a.y-fty-a-jS+i;!-*), 



(26) 



(27) 



where 



r(a) T(^) r(y+n) 

and definitions of these solutions on the whole complex plane are not shown here. 

They are connected by Kummer connection formulas in common converging 
regions, which include 

r(l-y)r(a + /S-y+l) r(y-l)r(a + /S-y + l) 

-3 = r(a _ y+1)r( ^_ y+1) - 1 + ^ (^) 



r(i-«)r(i- 



r(=<-7+i)r(/?- r +i) 



Logarithmic solution near z = 0 for y —> 0. To prove Eq.(13), we only need to 

Wi 

guarantee that its numerator — a^W2 approaches zero when y 0. 

r(yj 

According to Eq.(26), when y — > 0, 

Wl ~ 1+ r ( ,) g /fef i r(a+ " +1)rtf+ " +1) r(} ' +2) Z " 
1 r(«+i) r(/; + i) r( 7+ „+2) z (31) 

=l + r(y)a^z 2 Fi(a+l,^+l;2;z); 



and 

n« 2 <»a2J'i(«+l./!+l;2;2); 
Then the Eq.(13) can be understood. 



(32) 
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